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Abstract. The kHz QPOs observed in a neutron star low mass X-ray binary are likely to be 
produced in the innermost regions of accretion disk around the neutron star The rotational dynamics 
of the inner disk can be characterized by the presence of either sub-Keplerian or super-Keplerian 
accretion flow depending on the relative fastness of the neutron spin as compared to the Keplerian 
frequency at the inner disk radius. Within the magnetosphere-disk interaction model, the frequency 
difference between the two kHz QPOs observed in a given source can be estimated to be slightly 
higher than or nearly around the neutron star spin frequency if the neutron star is a slow rotator and 
less than the stellar spin frequency if the neutron star is a fast rotator. 

(~| . Keywords: Accretion and accretion disks. Oscillations, Neutron stars. X-ray binaries 

Pj PACS: 97.10.Gz,97.10.Sj,97.60.Jd,97.80.Jp 

INTRODUCTION 

The observations of neutron stars in low mass X-ray binaries (LMXBs) led to the 
discovery of the kHz quasi-periodic oscillations (QPOs) and the burst oscillations [8]. It 
has almost been firmly established from the observations that the frequency of the burst 
oscillation is very close to the spin frequency of the neutron star, i.e., Vburst — Vspin [5]. 

The kHz QPOs from neutron star LMXBs have some important observational proper- 
ties. These high-frequency QPOs usually appear in pairs in the c:^ 200 — 1200 Hz range. 
Both the upper kHz QPO frequency V2 and the lower kHz QPO frequency Vi correlate 
with the X-ray luminosity of the source [8]. According to early observations, the peak 
separation between the kHz QPO frequencies seemed to be directly related to the spin 
frequency of the neutron star, i.e., Av = V2 — Vi ~ Vburst — Vspin- The systematic analysis 
of all QPO data, however, showed that the peak separation Av is not constant for a given 
source. Indeed, Av decreases as both the kHz QPO frequencies V2 and Vi increase [6]. 
According to the analysis of the early kHz QPO data, Av ~ Vburst — Vspin for Vspin < 400 
Hz and Av ~ Vburst /2 — Vspin/2 for Vspin ^ 400 Hz among neutron star sources with dif- 
ferent burst or spin frequencies [9]. The recent analysis of all available data, however, 
revealed that Av significantly deviates from either Vspin or Vspin/2 for almost half of all 
sources and that A v/ Vspin decreases from a value > 1 to a value < 1/2 as Vspin increases 
[5]. 

The frequency correlations of kHz QPOs can be explained, albeit qualitatively at this 
stage, within the boundary region model [1, 4]. Next, we discuss the interpretation of 
the observational properties of kHz QPOs within the boundary region model. 



BOUNDARY REGION MODEL 



According to the boundary region model (BRM), the likely origin of the high-frequency 
QPOs is the innermost region of the accretion disk around the compact object. In 
neutron star sources, the kHz QPOs are produced in the boundary region (BR) or the 
boundary layer where the angular velocity of the accretion flow is non-Keplerian [1,4]. 
In black hole candidates, it is the relativistic disk region near the innermost stable 
circular orbit which is responsible for the production of the high-frequency QPO pairs 
with a frequency ratio around 1.5 [2]. 

The interaction of the inner disk with the magnetosphere of the neutron star rotating 
with an angular frequency Clspm = leads to the formation of a magnetohydrodynamic 
BR throughout which the angular frequency Q. of the accreting matter deviates from the 
Keplerian frequency to match the stellar spin frequency 12* at the innermost disk 
radius [3]. 

In the conventional regime of accretion where the neutron star is a slow rotator, 
Q.^ < QkIHii) and the angular frequency profile of the matter is characterized by a 
maximum frequency Q-mux = ^^('"o) at some radius ro in the BR [3]. In the Newtonian 
regime, the radial epicyclic frequency fc is the same as the angular frequency D. if 
Q. = Q.K- It is a well known fact that the nondegeneracy between k and CI appears 
in the regime of strong gravity due to general relativistic effects. Even in the non- 
relativistic regime, the degeneracy between K and Q is removed in a hydrodynamic 
BR throughout which fc > ^2 if the BR is sub-Keplerian [1, 4]. The analysis of global 
hydrodynamic modes of free oscillations in a typical sub-Keplerian BR revealed that 
the fastest growing modes are excited near the disk radius ro where Q. is maximum. 
The frequencies of the growing hydrodynamic modes match the test-particle frequency 
branches k±Q. and k in the limit of small hydrodynamic corrections. The same modes 
grow in the ~ rjn — ro range. The difference between the two consecutive frequency 
bands of the growing modes is ~ H which corresponds to the frequency separations 
Aco ~ f^max — 1.2— 1.3^2* and Aco ~ at r = ro and r = rjn, respectively [4]. Note that 
the estimations of the BRM for the peak separation of the growing mode frequencies 
are in agreement with the distribution of the Av/Vspin values observed for the relatively 
slowly rotating neutron stars with spin frequencies below ~ 400 Hz [5]. The frequencies 
of the modes which grow in amplitude at any particular radius in the BR increase all 
together as the local sound speed Cs increases [4]. The magnetic field strengths of neutron 
stars in LMXBs are generally thought to be weak enough for the radiation pressure to be 
important in the innermost disk regions. For a radiation pressure dominated inner disk, 
C.S o<: M, where M is the mass inflow rate in the inner disk [7]. We expect, in accordance 
with the correlations of the kHz QPO frequencies, that the growing mode frequencies 
correlate with the X-ray luminosity Lx °^ M. Figure 1 exhibits the correlation between 
the frequencies Vi and V2 of two consecutive modes growing at the innermost disk radius 
rin and the correlation between the frequency difference Av = V2 — Vi and the upper 
mode frequency V2 of the same modes estimated by the BRM [4] for a putative neutron 
star with a spin frequency 300 Hz. Note from Figure 1 that the power-law fit to the kHz 
QPO data of Sco X-1 [6] can also be used to fit the sample data estimated by the BRM 
simply by adjusting the proportionality constant in the power law. 

According to the present distribution of Av/Vspin over Vspm [5], the relatively fast 
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FIGURE 1. Power-law fits to the correlation of the frequencies Vi and V2 and to the correlation between 
the frequency difference Av and V2 estimated by the BRM for the growing modes at the innermost disk 
radius. The power law is given by Vi = C(V2/1000)^'^ Hz, where C = 795 for the solid curve and C = 724 
for the dashed curve [6]. 



rotating neutron stars with Vspin ^ 400 Hz seem to exhibit kHz QPO pairs for which Av 
is around or sometimes even less than Vspin/2. The systematic trend of observing small 
values of A v / Vgpin for sufficiently high values of Vgpin could be the result of a change 
in the structure of the BR. In the magnetosphere-disk interaction model, the structure 
of the BR depends on the fastness parameter, co* = ^2*/^2K(^n)• The neutron star is a 
slow rotator if CQ^ < 1 . For a slow rotator, the BR is sub-Keplerian as mentioned above. 
In the accretion regime where the neutron star is a fast rotator (co^ > 1), we expect to 
find a super- Keplerian BR for < r < ta, where ta is the effective magnetic coupling 
radius, e.g., the Alfven radius. The magnetic coupling radius, for a fast rotator, exceeds 
the corotation radius rco, where Q.^ = Q.}^. The angular frequency Q. of the disk matter 
accreting through the super-Keplerian BR deviates from the Keplerian frequency to 
match the stellar spin frequency Q.^ at r = rjn. Condition for the fast rotator to be an 
accreting neutron star and not to be a propeller can be written as Q.k <Cl< V^Q-k for 
fin^f^ fA- Thus, a fast rotator accretes only within a limited range of the fastness, i.e, 
I < CO^ < y/l. In Figure 2, we display the run of the dynamical frequencies K (dashed 
curve) and D. (dotted curve) throughout the super-Keplerian BR for a fast rotator with 
(0* = 1 .4. The angular frequency H lies in the — V2£Ik range (region bounded by 
the solid cuves) as expected. The radial epicyclic frequency K is always less than 
and sometimes even less than As compared to the right panel of Figure 2, r^ has 

a greater value in the left panel. Note that K'(rin) decreases as ta decreases. Although 
the mode analysis is necessary to distinguish among different hydrodynamic modes, it 
is noteworthy to find the frequency separation of two consecutive modes to be Aft) ~ K 
if the high-frequency modes with Q ± fc and Q. branches determine the kHz QPOs in the 
fast rotator regime. For ihtfast rotators. Aft) decreases as M increases since ta °^ M^'^l'' 
for the Alfven radius. 
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FIGURE 2. The radial profiles of the dynamical frequencies K (dashed curve) and CI (dotted curve) in 
the BR of the accretion disk around a rapidly rotating neutron star with a fastness (ft* = 1.4. The two solid 
curves correspond to the frequencies Q.k and V2Q.K- 



CONCLUSIONS 

The BRM estimations for the peak separations and correlations of kHz QPOs are in 
agreement with observations. A sub-Keplerian BR is appropriate for the slow rotators 
that exhibit kHz QPOs with the peak separations that are sUghtly higher than or nearly 
around or slightly less than the neutron star spin frequency. Within the magnetosphere- 
disk interaction model, super- Keplerian BRs are more appropriate for the accretion disks 
around fast rotators. In a super-Keplerian BR, the frequency difference of the mode 
branches is always less than the stellar spin frequency if the neutron star is a fast rotator. 
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